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ABSTRACT 


Thrust chamber performance should be evaluated in terms of an 
analytical model incorporating all the loss processes that occur in a real 
rocket motor. The Performance Standardization Working Group of JANNAF 
has identified the important loss processes in the real thrust chamber, and 
has developed a methodology and a recommended procedure for predicting 
real thrust chamber vacuum specific impulse. 

Simplified equations, based on the JANNAF reference procedure for 
calculating vacuum specific impulse, are developed to relate the delivered 
performance (both vacuum specific impulse and characteristic velocity) to 
the ideal performance as degraded by the losses corresponding to a specified 
list of loss processes. These simplified equations enable the various per- 
formance loss components, and the corresponding efficiencies, to be quanti- 
fied separately (except that interaction effects are necessarily arbitrarily 
assigned in the process). 

The loss and efficiency expressions presented can be used to evaluate 
experimentally measured thrust chamber performance, to direct develop- 
ment effort into the areas most likely to yield improvements in performance, 
and as a basis to predict performance of related thrust chamber 
configurations . 
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I. INTRODUCTION 


The JANNAF Performance Standardization Working Group started 
work in 1965, and in 1968 published the ICRPG Liquid Propellant Thrust 
Chamber Performance Evaluation Manual (Ref. 1) and several related 
reports. The manual presented a basic physical model for the rocket thrust 
chamber, including a listing of the primary loss processes, and described 
the methodology and computer programs which were available at that time 
for computing thrust chamber performance. 

Since 1968, the Working Group has continued to refine the methodology 
and the capabilities of the computer programs. In particular, an earlier 
empirical approximation used to describe the energy release loss is being 
replaced by a model based on physically describable processes. New com- 
puter programs have also been developed to implement the improvements in 
the methodology. 

New and revised manuals describing the improved performance pre- 
diction methodology and its application to correlation and analysis of mea- 
sured thrust chamber performance are in the process of preparation and 
adoption by the Performance Standardization Working Group, but will prob- 
ably not be published and distributed until some time in 1973. 

This report describes approximate procedures that can be used to 
correlate and evaluate experimentally measured thrust chamber perform- 
ance during the period between the general acceptance of the improvements 
in the methodology by those working in the field and the publication of for- 
mally ac 'ed manuals describing the improvements and their application. 
The appro imate procedures described here are completely consistent with 
the physical model and the loss process descriptions currently envisioned 
for the revised JANNAF reference methodology (see next section). 
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Two rocket motor performance parameters are considered. The 
vacuum specific impulse I S p vac is the thrust per unit mass flow of propel- 
lant. This parameter directly determines the performance and payload of a 
rocket-propelled vehicle. The vacuum specific impulse must be measured 
or predicted with the highest possible accuracy; the current goal is ±1% for 
new propellants and conditions, and even closer prediction for well- 
characterized systems. Both the experimental measurements and the analyt- 
ical techniques used in evaluating the vacuum specific impulse are capable of 
achieving this accuracy goal in most instances. 

The characteristic velocity c* is a mass flow parameter. It is used 
both to predict mass flow (which, in conjunction with the vacuum specific 
impulse, determines thrust level) and to obtain the preliminary estimation of 
the effect of the various loss processes on performance. (Note that the 
effect of a given loss process on characteristic velocity is not exactly the 
same as its effect on vacuum specific impulse. ) The characteristic velocity 
can be determined experimentally by testing a motor having a low exit-area- 
ratio exhaust nozzle; thrust measurement is not required. Experimental 
determination of c* is generally less accurate than I gp (about ±2 or 3% 
uncertainty) because of various measurement and data correction problems 
and assumptions involved in the definition of this parameter. An accuracy of 
±2 or 3% is considered satisfactory for design purposes, since mass flow 
rates can usually be adjusted by this amount to obtain a desired thrust level 
in a vehicle. Both the vacuum specific impulse and characteristic velocity 
are relatively insensitive to small changes in the level of the chamber pres- 
sure at which they are evaluated. 

Vacuum specific impulse is the preferred parameter for performance 
evaluation and correlation, because specific impulse is more meaningful and 
can be determined with better accuracy than the characteristic velocity. 

The analytical model employed for predicting thrust chamber perform- 
ance can be used to show the effect on the potentially achievable performance 
of each of the included loss processes. This capability is useful during 
rocket engine development, as it enables the magnitude and acceptability of 
individual loss components to be evaluated and thus shows where further 
development work would be profitable. 


2 


JPL Technical Memorandum 33-548 


As a first principle, experimental performance is to he stated as of 
the conditions under which it is measured and the reference ideal perform- 
ance calculated to correspond to the specified test conditions before compari- 
son with the test data. 

In using the approximate methods of performance prediction, all eight 
of the losses and effects listed in the next section must be considered and 
accounted for. (Note that any unknown or unaccounted-for losses or effects 
will be forced to appear as part of one of the recognized losses during the 
correlation process. ) In addition, it must be recognized that the approxi- 
mate methods result in an arbitrary distribution of the interaction effects 
among the various losses. 

Five primary steps are involved in the analytical correlation and 
evaluation of experimentally measured thrust chamber performance (see 
Fig. 1): 

(1) Determination of the experimental performance parameters for 
the actual operating conditions, based on specific measured data 
and on appropriate corrections to the measured data. 

(2) Specification of the conditions under which the experimental per- 
formance has been obtained. 

(3) Calculation of the analytically predicted performance correspond- 
ing to the specified test conditions, with the implicit determina- 
tion of the magnitude of all losses. 

(4) Comparison of experimental and predicted performances. The 
predicted and the experimentally measured performances are 
compared to see if they are in acceptable agreement. The agree 
ment criteria will depend upon the uncertainty associated with 
both the predicted and the measured performances; but agree- 
ment to within ±1% is the desired goal. If acceptable agreement 
is not achieved initially, the test measurements, input data, 
modeling, and calculations must be reviewed and errors found 
and corrected. 1 This process must be repeated until acceptable 


*The reference thrust chamber performance model can be changed only by 
action of the full JANNAF Performance Standardization Working Group. 
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agreement (correlation) is achieved. It is then presumed that 
the loss magnitude and efficiencies associated with the corre- 
lated performance prediction are a correct representation of the 
losses occurring in the experimental thrust chamber. 

(5) Evaluation of performance and losses. The overall perform- 
ance, the overall performance efficiency, and the individual 
losses and their corresponding efficiencies are evaluated to 
determine whether they are acceptable, or whether improvement 
in any area appears feasible with the state-of-the-art. At this 
point, a decision is made to continue development or to accept 
the thrust chamber design and performance at its existing stage 
of development. 

The correlation and evaluation procedure described above, and 
charted in Fig. 1, serves as confirmation of the satisfactory completion of a 
thrust chamber development program, or as a guide in determining the areas 
where further development work needs to be done. 

II. JANNAF PERFORMANCE PREDICTION METHODOLOGY 

The JANNAF performance prediction methodology currently enables 
calculation of steady-state vacuum specific impulse 

I 

sp 

vac 

pred 

of liquid-liquid propellant injection thrust chambers which (1) have only 
gaseous combustion products, (2) have conventional de-Laval nozzles, and 
(3) are large enough that the boundary layer influences only a small part of 
the total flow. The methodology takes into account the following interacting 
losses, with reference to one -dimensional isentropic equilibrium flow per- 
formance as ideal: 

(1) Energy-release loss, consisting of two parts: 

(a) Vaporization loss - due to incomplete liquid droplet vapori- 
zation at the nozzle throat, including thrust effects 
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due to both reacted gas and remaining liquid droplet 
components. 

(b) Mixture-ratio distribution loss — due to macroscopic non- 
uniform lateral distribution of the local time-averaged 

2 

mixture ratio of vaporized propellant at the nozzle throat. 

(2) Kinetic loss — due to rate-limited equilibrium shift during 
expansion. 

(3) Two-dimensional loss — due to upstream nozzle throat curvature 
(which affects only mass flow and, hence, c*), and to nozzle 
divergence shape and exit angle (which affect only the vacuum 
specific impulse). 

(4) Boundary-layer loss — due to friction and heat transfer at the 
wall. 

Within the context of, or as input to, the computer programs used in the 
methodology, it is possible to account for effects due to 

(5) Propellant impurities. (Adjust enthalpy and composition of 
injected propellant. ) 

(6) Heat lost to the thrust chamber from the injector face and the 
chamber wall upstream of the boundary -layer attachment point. 
(Reduce enthalpy of injected propellants. ) 


^There can be an additional mixture-ratio distribution loss due to timewise 
variations from the mean in the mixture ratio of the vaporized propellent 
flowing through any element of the nozzle throat cross section. This 
temporal component of the mixture-ratio distribution loss is not presently 
accounted for in the JANNAF methodology. 

3 In some rocket motors, there is a turbulent region just downstream from 
the injector, followed by a region in which the boundary layer appears to 
develop normally. Heat transfer upstream of the effective starting point for 
the boundary layer and radiant heat transfer throughout the thrust chamber 
must be accounted for separately from the boundary -layer calculations. 

The location of the elective starting point for the boundary layer is deter- 
mined experimentally by examination of the axial variation of the wall heat 
flux. 
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( 7 ) 


Mass addition into boundary layer. (Use MABL or BLIMP 

,4 

programs. ) 

(8) Feed-system energy loss and exchange, such as that due to 
pumps, turbines, and regenerative-cooling heat transfer. 

(Adjust enthalpy of injected propellants. ) 

In order to develop the analytical expressions for predicting the 
delivered performance, the real rocket thrust chamber is modeled as a num 
ber of inviscid stream tubes flowing within a virtual nozzle. The virtual 
nozzle is offset from the real nozzle wall by a distance 6"'', obtained from 
the boundary -laye r analysis, such that the mass flux of the real and the 
inviscid flows is identical (Fig. 2). There is no mixing across stream-tube 
boundaries . 

The stream tubes are constrained to flow together through the throat 
of the virtual nozzle. The individual and the total stream-tube mass flow 
are related to the throat size and curvature, and to the pressure level in the 
combustion chamber. 

The thrust chamber vacuum specific impulse 


sp 

vac 

pred 


is obtained by summing the thrust of the equivalent inviscid-flow stream 
tubes and adding corrections for the pressure force and the momentum 
deficiency of the flow in the boundary layer at the nozzle exit, and dividing 
by the discharged mass flow rate. 

For liquid-liquid spray injection, the composition of the gases flowing 
through the stream tubes is given by a distributed energy -release combus- 
tion model (DER), which assumes that droplet distribution from the injection 
elements establishes the stream-tube configuration and that droplet evapora- 
tion is the controlling factor in the combustion energy-release process. 

Thus, a stream tube contains at any location a quantity of completely reacted 


*Mas s addition boundary-layer programs being developed by the JANNAF 
Performance Standardization Working Group. 
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gas and a residual of unevaporated liquid droplets; the mixture ratio of the 
evaporated and reacted combustion gas may differ from the initial injected 
mixture ratio of the stream tube. 

The basic JANNAF performance prediction methodology computes 
vacuum specific impulse by the following procedure: 

(1) The distributed energy-release series of computer programs is 
used to compute the vaporized (and reacted) mass, mixture- 
ratio, enthalpy, and stagnation pressure (or entropy) distribu- 
tions of the flow at the nozzle throat. Input to this program is 
basically the injected propellant mass, enthalpy, mixture -ratio, 
and droplet-size distributions. 

(2) The TDK program uses the output of the DER program to compute 

the two-dimensional kinetic multi-stream-tube inviscid-flow noz 

zle performance. This initial calculation, TDK1 , is made with 

the geometric nozzle dimensions, in order to obtain boundary 

stream-tube properties that can be used to compute boundary- 

5 

layer parameters. 

(3) A boundary -layer program, either an integral-method program 
such as TBL or a finite-difference-type program such as MABL 
or BLIMP, is used to compute the boundary-layer momentum 
and displacement thickness, using input data from the TDK 
program. 

(4) A final TDK calculation, TDK2, is made using the virtual nozzle 
dimensions, to obtain the inviscid gas -flow thrust and the prop- 
erties of the boundary stream -tube flow at the nozzle exit. 


5 Note that TDK is a package of programs which also contains the one- 
dimensional equilibrium program ODE, the one -dimensional kinetic pro- 
gram ODK, a transonic program TRANS, and a TDE option. A two- 
dimensional equilibrium program TDE may be added to the programs 
available, either separately or as an option in the TDK program package. 
TDE may be useful as a less -expensive substitute for TDK in obtaining 
boundary stream-tube properties for computing boundary -layer parameters. 
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(5) The total thrust is the sum of 

(a) The thrust resulting from the inviscid kinetic expansion of 
the throat-plane vaporized and reacted gases through the 
effective supersonic nozzle contour (from TDK2). 

(b) The momentum flux of the unevaporated liquid droplets 
passing through the nozzle throat plane (from DER). 

(c) The thrust deficiency of the real gas flow, with its bound- 
ary layer, compared to the assumed inviscid flow (from 
TBL, etc. , and TDK2). 

At present the assumption is made that interaction of the already- 
combusted gases with the remaining liquid droplets (momentum exchange, 
evaporation, energy release) downstream of the throat plane has no net 
effect on the thrust as given by items (5)(a), (b), and (c), above. The per- 
formance losses and the computer programs used in calculating performance 
are summarized in Tables 1 and 2. 

Modifications and extensions to the methodology to enable prediction of 
performance with gas -liquid and gas -gas propellant injection are currently 
being considered, and at some time in the future it is anticipated that pro- 
visions will be made to handle 2 -phase combustion products. 

It is important to recognize that, at the present stage of development 
of the JANNAF performance evaluation procedures, only loss mechanisms 
(1) through (8), above, are considered. Meaningful correlation and evalua- 
tion of experimental data cat. be achieved only when the comparison analyt- 
ical model incorporates all of the los3 processes present in the real motor 
being evaluated. Fortunately, the losses (1 through 8) inventoried are ade- 
quate to model most currently important thrust chamber configurations. 
Thus, the JANNAF methodology serves a useful purpose in its present state 

of development, ven while steps are being taken to extend it to cover more 

6 

complicated systems. 

6 The computer programs ODE, ODK, TDK, TBL, and MABL can be 
obtained from CPIA, 3261 Georgia Ave. , Silver Spring, Maryland 20910 
(Attn: Tom Reedy). The computer program DER, in a version which 
handles both liquid-liquid and gas -liquid injection, has been developed by 
Rocketdyne for JPL under Contract NAS 7-746, and is currently being 
checked by JPL. In its current form, DER does not appear to model 
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Ill EXPERIMENTAL MEASUREMENTS AND SPECIFICATION 
OF OPERATING CONDITIONS 


If experimentally measured and predicted thrust chamber performances 
are to be correlated, it is necessary that all of the required measurements 
be made to enable determination of the experimental performance and that 
the conditions under which the performance is measured be completely speci 
fied, so that the corresponding performance prediction can be made. 

The experimental data, which are required to determine delivered 
values of vacuum specific impulse and characteristic velocity or as input to 
the performance prediction procedure and which must be measured during 
thrust chamber operation, are listed below and indicated in Fig. 3. 


^amb 

measured thrust 

m o’ m f 

propellant flow rates, oxidizer and fuel 

T. P 

propellant temperature and pressure, oxidizer and 
fuel, at injector inlet 

A t’ A e’ A lip 

nozzle throat, exit, and lip areas (Since it is difficult 
to measure these during motor operation, corrected 
"cold" measurements are usually used. ) 


static pressure at a specified axial location and 
chamber area ratio 


ambient pressure 

p lip 

"base pressure" at nozzle lip 

T , q 

w w 

temperature and heat flux profiles at chamber and 
nozzle wall (Measured values are desirable, though 
these quantities can be calculated. ) 


adequately the gas -liquid injection case, and improvement will have to 

await results of fundamental studies in this area. t V'^TmTC CPIA bv 
existing DER program will become available through COSMIC or CPIA by 
March 1973. Improved versions of DER will print out interface informa- 
tion needed to start the TDK calculation, and both DER and TDK will print 
out information needed to start boundary -layer calculations with TBL, 
MABL or BLIMP. BLIMP is an alternate boundary -layer program which 
is currently being modified to fit into the JANNAF family of programs. 
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Additional information, beyond the data which must be obtained during 
testing, is needed as input for the performance prediction computer pro 
grams and procedures, and for the calculation of the experimental delivered 
performance. This information fall, into three general categories: 

(1) Chamber and nozzle size and shape, including: 

(a) Chamber diameter and length. 

(b) Location and area ratio at chamber static pressure taps. 

(c) Throat radius of curvature. 

(d) Nozzle throat area, exit area, and lip area. ("Cold" 
measurements should be corrected for thermal and pres 
sure effects to obtain effective "operational" values. ) 

(Z) Propellant inlet conditions : 

(a) Composition, including impurities. 

(b) Enthalpy, including any effects due to regenerative cooling 
and impurities. 

(3) Propellant injection characteristics: mass, mixture ratio, drop 

size, and enthalpy distribution. 

IV. EXPERIMENTAL DETERMINATION OF VACUUM 
SPECIFIC IMPULSE 

Vacuum specific impulse can be determined directly from tests made 
in a vacuum environment, or it can be calculated from measurements made 
during test of a low-area-ratio configuration (fully attached nozzle flow) 
thrust chamber in an ambient pressure environment. The physical and geo- 
metrical data, which are required both to determine the experimental value 
of the vacuum specific impulse and to define the corresponding analytical 
model, are listed in the preceding section. 
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The general expression for the experimentally measured thrust 
chamber vacuum specific impulse is 


I 

sp 

vac 

exp 


F , + p A + /(p _ p, • ) dA 1 . 

amb K a e J ^a Hip lip 


m. 


( 1 ) 


This expression includes corrections for the external pressure effect on the 

lip thickness of the real nozzle, a factor which, while usually small, should 

7 

not be neglected. The quantities are defined in Fig. 3. 


At this point, it is convenient to define an overall thrust chamber 
specific impulse efficiency 


TC 

exp 


I 

sp 

vac 

exp 


I 


sp 

vac 


ODE 

(0/F)avg 


( 2 ) 


^ An overall specific impulse for an entire engine assembly can be defined as 


SF 

*sp ~ E(mass discharged) 

engine 

where (1) EF includes all thrust contributions, from the main thrust 

chamber and from other sources, such as turbine exhausts 
and separately exhausted coolant flows. 

(2) x(mass discharged) includes the mass discharged from the 
main thrust chamber, that discharged from auxiliary pro- 
pulsive dumps, and that which is vented or leaks from any 
part of the system whether or not it produces any net 
thrust. 


This report deals only with the thrust chamber psrformance parameters, 
and does not consider "external ' effects associated with feed systems, 
pumps, turbine exhausts, vents, or thrust chamber coolant flows which are 
expanded and discharged separately from the combustion chamber gases. 
Any references to system" within this report refer only to the thrust 
chamber and its regenerative cooling flows, if any. 
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where 


I 

sp 

vac 

ODE 

(0/F)avg 


is the ideal thrust chamber performance that would be attained in the absence 
of all of the real thrust chamber internal loss processes. It is calculated 
using the ODE computer program for the actual thrust chamber test 
injector-inlet propellant composition and enthalpy, and nozzle e::it-area ratio 
Note that, for a regene ratively cooled chamber, the injector-inlet propellant 
enthalpy must include the heat transferred to the propellant from the 
chamber, and is not the propellant enthalpy at the thrust chamber inlet. 

V SIMPLIFIED FORM FOR ANALYTICAL PREDICTION OF 
VACUUM SPECIFIC IMPULSE 


A close analytical approximation to the JANNAF reference procedure 
for calculating thrust chamber specific impulse is 


sp 

vac 

pred 


sp 


1 

vac 

TDK u 

(0/F)f 
' 'i, vap 

HL 

€ ' 

drop interact 


| m i, vap \ . A 1 F BL 


AF 


drops 


(3) 


rrr 


where 
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/ 1 \ and (0/F)f are stream-tube evaporated mass and 

' ’ i.vap i.vap 

mixture ratio at the throat. 

(2) A F is a thrust decrement which accounts for friction and 
heat-transfer effects in the boundary layer, based on flow through 
the equivalent inviscid-flow nozzle of area ratio t . 

(3) AF is a thrust increment equal to the momentum flux of the 

liquid droplets passing through the throat plane of the nozzle. 
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(4) The subscript HL indicates that the enthalpy of the propellant at 
the injector inlet has been reduced to account for heat loss 

in addition to that accounted for by the boundary-layer computer 
program (TBL, MABL, or BLIMP). The additional heat losses 
for which correction should be made are 

(a) Heat lost to the system by convective or radiative heat 

transfer in the region upstream of the turbulent boundary- 
layer start point. This heat can be stored in the thrust 
chamber wall, lost to the outside environment by convec- 
tion or radiation, or picked up by a regenerative coolant 
in such a manner that the measured propellant enthalpy at 
the injector inlet is increased. (Heat which is picked up by 
the propellant but does not increase measured enthalpy at 
the injector inlet, for instance heat recycled from the 
injector face plate, does not enter into this correction. ) 

(b) Heat losses due to radiation from the combustion gases to 
the chamber and nozzle walls in the turbulent boundary 
layer region, and, in the divergent portion of the nozzle, 
the heat lost by direct radiation from the gases to the out- 
side environment. 

(5) The subscript €' indicates that the quantities are determined 
for an equivalent inviscid-flow nozzle, obtained by displacing the 
boundary of the real nozzle by the displacement thickness of the 
boundary layer 6* . 

(6) The subscript "drop interact" indicates that the enthalpy of the 
vaporized propellant has been reduced to account for the energy 
used to accelerate the liquid droplets remaining at the nozzle 
throat. 

The radiative heat loss from the nozzle-exit section can be significant 
for some small thrust chambers, but no standard method of accounting for 
this loss has been developed, and it is not considered further in this paper. 

Heat radiated to the inside walls of the thrust chamber can be treated 
as an addition to the convective heat transfer (as far as total heat loss and 
heat balance are concerned) in both the region upstream of the turbulent 
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boundary-layer start point and in the turbulent boundary-layer region itself, 
but will not enter into the calculation of the boundary-layer displacement and 
momentum thicknesses other than as it increases the temperature of the 
inside surface of the combustion chamber. 

The various heat losses and their corresponding treatment in the per- 
formance prediction procedure are indicated in Fig. 4. 

The droplet thrust term will be ignored from here on in the develop- 
ment and discussion of the simplified performance prediction procedure. 

This can be justified by the following arguments: 

(1) In any acceptable rocket motor, the unevaporated propellant drop- 
let fraction will be small, and the major part of the residual 
droplet effect will show up in the decreased mass of the reacted 
product gases . 

(2) The thrust of the droplets at the throat plane arises from 
entrainment in and momentum transfer with the surrounding 
reacted gases, with a corresponding slowing down of the gases. 
For small droplet fractions, it is reasonable to offset the throat 
momentum flux of the droplets against the associated slowing 
down of the gases at the throat. This is effectively accomplished 
by using the existing computer programs (ODE, ODK, and TDK) 
with the as -injected propellant enthalpy to calculate the reacted 
gas properties and performance, and omitting the droplet thrust 
term. 

Equation (3) can now be replaced by Eq. (4), where the stream-tube 
performance is calculated for single-phase gas flow through the equivalent 
inviscid-flow nozzle, the droplet thrust term has been omitted, and the 
boundary-layer term has been expanded in terms of the boundary -laye r 
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parameters 6 and 6* and the gas properties in the boundary-layer stream 
0 0 ^ 
tube at the exit station of the equivalent inviscid -flow nozzle: 


X sp ».p. 

pred 


?%“)• 


iZp A cos o 
1 e e 






(4) 


TDK 

€' 


i 

TDK 

(O/F)* 

i,vap 

HL 
€ ' 

Additionally, the subscript "vac" has been dropped, and it is from here on 
to be understood that specific impulse is always calculated for vacuum 
conditions . 

It is desirable that a simplified performance prediction calculation 
procedure be referenced to the real nozzle exit-area ratio and dimensions. 
The exit-area ratio of the real and of the equivalent inviscid-flow nozzles 
are related by the expression 
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Now Eq. (4) can be pul in the form 
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8 The boundary-layer parameters 0 e and 6* are defined and discussed in 
Appendix A. 
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From this point on all quantities will be evaluated at the geometrical nozzle 

exit-area ratio t. The fluid properties in the last two terms of Eq. (6), 

(o u and p ) are understood to be evaluated at the edge of the equivalent 
Nr e’ e e 

inviscid flow at the exit plane of the nozzle (TDK). 

Since the second term of Eq. (6) is a change in thrust associated with 
a change in nozzle exit-area ratio caused by the boundary layer, this term 
is combined with the third and fourth terms as the total boundary-layer 
effect on specific impulse, ^2 I sp BL . to be used when specific impulse is 
referenced to the value corresponding to inviscid flow through the geo- 
metrical nozzle. 

The predicted specific impulse is then given by 
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The next step is to develop simplified approximations ana methods of 

evaluation for both the first and second terms of Eq. (7). The evaluation of 

the boundary -layer loss, A£lsp , is treated separately in Appendix A, 

BL. ^ 

As a first step, Eq. (7) can be written as 
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9 A somewhat similar development to that following, but using subtractive 
efficiency increments, is described in Ref. 2. 
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and 
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The chamber injection-end heat loss factor 


sp 

HL 


accounts for the convective heat lost to the system from the region upstream 
of the starting point for the boundary layer (injector faceplate and a portion 
of the chamber wall) and for radiation loss throughout the chamber; this is 
the heat lost in addition to that accounted for by the boundary-layer calcula- 
tions. This factor is less than unity only for that portion of the heat which is 
lost to the system or transferred to a regenerative coolant so that the tnea- 
I^Ted injection enthalpy is increased; there is no net performance loss when 
the heat is picked up internally by the propellant (as in an injector manifold 
and faceplate) and returned to the combustion region. The treatment of the 
various heat losses as they affect performance is shown in Fig. 4 The heat 
loss is converted to an enthalpy loss by dividing by the total propellant flow 
rate. The sensitivity of the specific impulse to enthalpy change can be deter- 
mined by computing ODE performance using perturbed values of the standard 
px'opellant enthalpies; the results can be plotted in the form of 
[( 1 / I )(AI /AH)] versus propellant mixture ratio for various nozzle exit 
area ratios, as shown in Appendix B, Fig. B-4. This factor can then be 
used, as shown by Eq. (10), to calculate the corresponding specific impulse 
correction factor. 
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The two-dimensional nature of the nozzle flow affects both the 
pressure at the nozzle exit and the net axial component of the exit momentum. 
It also affects the mass flow rate through the nozzle throat and the distribu- 
tion of pressure and pressure -dependent properties in the boundary-layer 
flow in the nozzle. The effect of two-dimensional flow on specific impulse 
is given by the factor 


2D 

This factor is given for several common nozzle configurations in Ref. 1, and 
results for conical nczzles are reproduced here as Figs. 5, 6, and 7. 

Equation (9) represents the first useful stage in the simplification of 
the performance calculation procedure; it substitutes the ODK computer pro- 
gram for the more complicated and expensive -to -use TDK computer program, 
by introducing the factor 


? D 

and it accounts for the heat loss in the chamber upstream of the boundary- 
layer attachment point by means of the factor 

11 1 

sp 

HL 

These factors can be completely or partially precomputed, usually using 
simplified approximations to the formal defining equations (Eqs. 10 and 11), 
and charts can be prepared for use in the simplified performance prediction 
procedures . 
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A further simplification of Eq. (9) is obtained by replacing each 
calculated ODK. stream-tube performance with an ODE. performance multi- 
plied by an appropriate kinetic efficiency factor. Thus, 
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In Eq. (12) each of the factors inside the summation sign is evaluated 
for the evaporated mixture ratio of the particular stream tube, and for the 
geometrical nozzle area ratio. 

The kinetic loss factor 



kin 


arises from the fact that the combustion gas composition cannot change fast 
enough to follow the equilibrium composition corresponding to the changes in 
pressure and temperature as the gases expand through the nozzle. This 
factor must be calculated for a specific propellant, operating conditions, and 
thrust chamber size and shape, as shown by Eq. (13), using the ODE and 
ODK computer programs. For a given pro ~ t and nozzle, values of 
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n I 

s ?i 

kin 

can be precomputed, as shown in Appendix B, Fig. B-3. 

Equation (12) represents the preferred simplified performance predic- 
tion procedure. It is the method that should be used if the stream-tube 
evaporated mass fraction and mixture -ratio distribution is given (as from a 
DER computer program calculation), or is provided from some other source, 
or can be reasonably estimated. 

The final stage in the development of a simplified performance equa- 
tion from Eq. (7) is to express the performance in terms of a series of 
efficiency factors, each of which represents the effect (mass -averaged over 
all stream tubes) of a specified physical loss process on the overall per- 
formance, with the performance referenced to the ODE vacuum specific 
impulse at the overall average mixture ratio. To do this it i3 necessary to 
introduce two new efficiency factors which account for the effects of propel 
lant vaporization and of mixture-ratio distribution on the overall perform- 
ance; the product of these two efficiencies is defined as the energy-release 
efficiency. 

The final form of the simplified performance equation then becomes 
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or the equivalent 


I 

sp 

pred 


n i ni nj n x 

sp sp sp sp sp 
HL 2D kin vap (0/F)dist 


-An 


i 

sp 

BL 


I 

sp 

ODE 

(0/F)avg 


€ 


TC 

pred 


I 

sp 

ODE 

(O/ F)avg 
€ 


(15) 


where 


sp 

TC 

pred 



sp 

HL 


n I 


sp 

2D 


sp sp 
kin vap 


sp 

(0/F)dist 


An x 


sp 

BL 


(16) 


In Eqs. (14), (15), and (16), 
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is a mean kinetic loss efficiency defined by the following equation: 
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The vaporization loss is defined by 
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and the mixture-ratio distribution, or stream-tube loss, is defined by 
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The two above-defined losses, which concern the uniformity of mixing 
and the completeness of droplet evaporation, are closely related and can be 
grouped together as the energy-release loss. Thus, 
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The vaporization loss factor (Eq. 18) accounts for the fact that only that 
fraction of the injected liquid which evaporates can react and release its 
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chemical energy to produce thrust. It also accounts for the fact that 
nonuniform droplet evaporation can cause vaporized mixture “ratio distribu 
tions (at the nozzle throat) which differ from the injected mixture -ratio 
distributions. 

The vaporized mass and mixture-ratio distribution needed to evaluate 



vap 


are obtained from the DER computer program calculation, which in turn 
requires as input a complete description of the injected mass, mixture-ratio, 
and droplet-size distributions. 

The primary mixture-ratio distribution loss (Eq. 19) arises from non- 
uniform injected mixture-ratio distribution. With the usual concave- 
downward variation of specific impulse with propellant mixture ratio, the 
mass -averaged performance of a distributed mixture-ratio flow is typically 
less than the performance which would be achieved with the entire flow at its 
average mixture ratio. 

The boundary -layer loss efficiency factor 
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used in Eqs. (15) and (16) is defined as 
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is defined by Eq. (8). The evaluation of the boundary -layer loss is discussed 
in more detail in Appendix A. 

The reference ideal performance 

I 

sp 

ODE 

(0/F)avg 

used in Eqs. (13) through (21) is evaluated for the specified test conditions of 
propellant composition (including impurities), injector inlet enthalpy, over- 
all mixture ratio, chamber pressure, and nozzle exit-area ratio. 

The ODE (one -dimensional isentropic equilibrium) performance is often 
available from calculations made for the propellants at standard tabulated 
initial enthalpy values. Performance at actual inlet enthalpy conditions can 
be obtained by correcting this standard enthalpy performance, using the spe- 
cific impulse -enthalpy influence factor developed to correct for heat loss 
from the combustion chamber walls (cf Eq. 10 and Appendix B, Fig. B-4). 

If the efficiency factors in Eq. (15) are evaluated according to the 
primary defining equations (Eqs. 17 through 23), Eq. (15) will yield the same 
predicted performance that would be given by Eq. (7) or Eq. (12). However, 
it must be recognized that, when the flow mixture -ratio distribution spans a 
region in which the performance (or an efficiency factor) is nonlinear with 
mixture ratio, apparent, but compensating, distortions will be introduced 
into related terms of Eq. (15). This is caused by using a reference specific 
impulse which corresponds to the average mixture ratio of a striated flow. 

Thus, if two stream tubes have mixture ratios below and above the mix- 
ture ratio for maximum performance, the reference specific impulse will be 
higher than the actual specific impulse of either stream tube, and the specific 
impulse efficiency for mixture-ratio distribution will have to be correspond- 
ingly low. However, if the stream-tube mixture ratios spanned a region in 
which the performance was nearly linear with mixture ratio, the specific 
impulse at the average mixture would be a ''real" value, and the specific 
impulse efficiency for mixture-ratio distribution would be near 100%. 

Similarly, an error can be introduced by choosing a specific impulse 
efficiency for kinetic effects at the average mixture ratio, rather than using 
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a mass -averaged value for the actual stream mixture ratios. Other 
distortions could be caused by the effect of vaporization efficiency in chang- 
ing effective stream -tube mixture ratios. 

Additional distortions in the relative values of the loss efficiency fac- 
tors, but no net error in the calculated performance, can result from the 
fact that there are interactions between the loss processes which cannot be 
correctly assigned when the losses are individually defined, as they must be 
for the simplified procedures developed herein. These effects cancel out 
internally when the efficiencies are defined according to the primary defini- 
tions given here. However, errors can be introduced when the approximate 
expressions for evaluating some of the loss efficiencies are used; probably 
the most important interaction involves energy release (vaporization) and 
kinetics. 

The above examples serve to indicate the necessity for care in using 
the form of the performance prediction equation given by Eq. (15). It is 
necessary to select values of specific impulse efficiency for kinetics, vapori- 
zation, and mixture -ratio distribution effects which correspond to the real 
stream-tube mass and mixture-ratio and vaporization distributions, or at 
least to reasonable estimates of these distributions. It should also be kept 
in mind that, because of the above-mentioned distortions, maximum pre- 
dicted performance may not correspond to maximum overall specific 
impulse efficiency. Thus, when two rocket motor systems are being com- 
pared, comparison should be on the basis of predicted specific impulse, 
rather than on overall efficiency. 

The above considerations emphasize again that simplified performance 
prediction by the procedure of Eq. (1Z) is preferable to that based on Eq. (15), 
even if it requires estimation or assumption of reasonable mass and mixture- 
ratio distribution and component vaporization efficiencies. 

Equation (15) can be used to make initial estimates of the performance 
of proposed thrust chambers, or to make parametric design studies around 
a given configuration. The two-dimensional-flow loss efficiency can be 
determined from existing design charts. The kinetic, boundary-layer, and 
upstream chamber heat-loss efficiency terms can be calculated by using the 
ODK, ODE, and TBL computer programs or charts prepared from these pro- 
grams. The energy-release efficiency can be calculated from postulated 
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injection conditions, using the DER computer program and Eq. (14). As a 
last resort, some or all of these quantities can be estimated on the basis of 
previous experience. 

If it is known that a given propellant system and a range of thrust 
chamber configurations are to be subjected to thorough analysis, it is worth- 
while to use the ODE, ODK, TDE or TDK, and boundary layer computer pro- 
grams early to investigate a broad parametric range of operating conditions, 
if feasible, the energy-release efficiency can be precalculated in the same 
way, using the DER and ODE computer programs. The results of these cal- 
culations can be used to develop a set of "influence coefficients, " which 
become the basis of a procedure for correcting design condition performance 

predictions to obtain predicted performance at actual test conditions or for 

r ... 10 
slightly modified design conditions. 

The results of a typical set of parametric calculations for the fluorine - 
hydrazine propellant in a given thrust chamber configuration are presented 
in Appendix B. 

Table 3 of this report summarizes the derivation of the simplified per- 
formance prediction equation and the definitions of the specific impulse loss 
process efficiency factors developed here. 

VI. CORRELATION OF EXPERIMENTALLY MEASURED AND 
PREDICTED VACUUM SPECIFIC IMPULSE 

The measured and predicted vacuum specific impulses can be com- 
pared to confirm or infer the magnitudes of the individual losses. The losses 
can then be evaluated and a judgment made as to whether they are individually 
reasonable and acceptable, or whether thrust chamber or injector design 
changes should be made in an attempt to decrease the magnitude of those 


10 These influence coefficients are sometimes used inversely to "correct 1 
experimental data to a common reference condition. Such a procedure 
violates the premise of this report that experiment and analysis should 
be kept completely separate and independent right up to the point of com- 
parison. However, it is convenient for removing secondary effects from 
a mass of experimental data so that primary effects can be more clearly 
recognized, and for facilitating comparison with predicted performance 
calculated at a "design" reference condition. 
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losses which are susceptible to control. Some losses, such as the 
two-dimensional divergence loss, the kinetic loss, and often the boundary- 
layer loss, are not very susceptible to the control. On the other hand, the 
energy-release losses (mixture -ratio distribution and vaporization) can be 
controlled, though usually at the expense of some other deliberately 
designed-for characteristic, such as low heat flux, chamber wall compati- 
bility, low injector pressure drop, or small combustion chamber size. 

If complete operating condition data are available, and it has been pos 
sible to use the DER, TDK, ODK, ODE, and boundary layer computer pro- 
grams to predict all of the loss components, then performance correlation 
consists of demonstrating that 
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where A is a correlation parameter defining an acceptable limit of 
cor 

error, based on the uncertainty involved in arriving at the two values ox 
specific impulse; a value of A cor = 0. 01 or better is suggested as a goal, 
though higher values may have to be accepted at times. 

If Eq. (22) is not satisfied at the first attempt, agreement within a 
reasonable limit of error can usually be obtained by recalculation after a 
careful examination for possible errors in experimental data, input informa- 
tion, and calculation procedures. Correlation in this manner is accepted as 
confirming that the magnitudes of the separate loss effects shown in Eq. (12) 
or (15) have been correctly predicted. 


11 


A more refined statistical approach to the correlation of measured and 
predicted performance which considers the uncertainty of each quantity 
will be given in the "JANNAF Performance Data Analysis Manual being 
prepared by Rocketdyne under Contract NAS8-28603. The basis for this 
approach is contained in Statistics , by W. L. Hays, published by Holt, 
Reinhardt, & Winston, 1963. 
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The same procedure is followed if some of the loss components have 
been determined from charts or even estimated on the basis of previous 
experience. The result in either case is a quantification of losses and 
corresponding efficiencies which can be used as a basis for evaluation of the 
thrust chamber design and performance. 

Figure 1 shows graphically the above correlation process and the sub- 
sequent procedure for evaluating the performance and the losses to deter- 
mine whether they are acceptable, or whether additional design and 
development work is needed. At a lower level of sophistication, the injection 
mixture -ratio and droplet-size distribution input data may not be known, so 
the DER computer program cannot be used to determine the energy-release 
efficiency components. In this case, the lumped energy-release losses can 
be "backed out" of the data (at a reduced confidence level of accuracy, 
because there is no overall check on the consistency of the results) by the 
following approximate formula, derived from Eqs. (16), (20), and (22), with 
the assumption that & CQr = 
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exp pred 
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In the above expression, it is presumed that the component loss effici- 
encies, other than 


ER 

can be calculated from the available data, and that actual performance has 
been measured. 

This procedure enables values to be assigned to all of the losses, so 
that the evaluation process shown in Fig. 1 can be completed. An example 
of a correlation and evaluation of test data obtained using the fluorine - 
hydrazine propellant in a given thrust chamber is presented in Appendix C. 
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VII PERFORMANCE PREDICTION FOR RELATED THRUST 
CHAMBER CONFIGURATIONS 


The individual loss components determined by correlating the test data 
obtained with a given thrust chamber can be used as the basis for predicting 
the performance of related thrust chamber configurations, using the same 
basic performance prediction and correlation methodology and equations. 

Related thrust chamber configurations may range from variations in 
propellant inlet temperature or changes in chamber wall cooling method and 
temperature to the addition of a high-area-ratio nozzle extension to a low - 
area-ratio test motor. Depending on the differences between the modified 
thrust chamber and/ or extrapolated operating conditions, it may be possible 
to retain a few of the correlated performance efficiency factors without 
change. The other performance efficiency factors, and perhaps the refer- 
ence ODE performance, mu3t be adjusted or recalculated according to the 
defining equations. Then the modified or extrapolated performance is calcu- 
lated from Eq. (12) or (15). 

It will be found that almost every case of performance extrapolation 
becomes a new problem in performance prediction, requiring recalculation 
of the reference ODE performance and adjustment of most of the efficiency 
factors. The recalculation task can be simplified if a sufficiently wide para- 
metric range of calculations is made initially for the propellant system under 
consideration, as described previously and shown in Appendix B. 

Appendix C illustrates the use of the calculated factors of Appendix B in 
correlating test data and in 'extrapolating" to performance at a different 
nozzle exit-area ratio. 

VIII. EXPERIMENTAL DETERMINATION OF 
CHARACTERISTIC VELOCITY 

The objective of this and the following sections of this report is to 
define the performance parameter, characteristic velocity c*, and to develop 
a procedure for predicting and correlating this parameter which is consistent 
with the JANNAF thrust chamber model and specific impulse prediction 
methodology. 
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The characteristic velocity is a mass-flow performance parameter. 

It is related to the energy release in the combustion chamber up to the nozzle 

throat, where the mass flow is determined. The definition of the experi- 

, . . 12 

mentally measured c* is 
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12 As defined, c* has units of [F t/M]. Since force is a derived unit which 
is defined in terms of the basic [L, M, t] units, c^\ as well as I S p» can 
be expressed in the alternate, but equivalent, [L</t| units. This is accom 
plished by multiplying by the unity unit conversion ratio corresponding to 
the measurement system being used. In the SI system this unity ratio is 


1 . 0 kg-m 

„ T 2 
N-sec 


In the English technical system of units, this unity ratio is 


32. 174 lbm-ft 
lbf-sec 2 

J 


Equivalent quantities of specific impulse and characteristic velocity in 
both the SI and the English technical systems are summarised in the 
following "conversion box": 
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where 


is a mean effective stagnation pressure at the nozzle throat. It 
is determined from the relation p* = p g ^Pq^s ^ 



is estimated from the well-known analysis for combustion in a 
cylindrical combustion chamber, depending on the tap location 
(see Fig. 8). 


A is the nozzle throat area during thrust chamber operation . It 

is usually determined by correcting the "cold" throat area for 
effects due to pressure, thermal expansion and thermal stress, 
erosion, etc. 


is a static pressure measured at the wall of the combustion 
chamber nozzle at the longitudinal station x. 


At this point, it is convenient to define an overall thrust chamber 
characteristic velocity efficiency: 
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Here 


c 5de 

(0/F)avg 


is the ideal characteristic velocity that would be attained in the absence of 
all of the real thrust chamber internal loss processes. It is calculated using 
the ODE computer program for the actual thrust-chamber -test injector -inlet 
propellant composition and enthalpy. Note that for a regeneratively cooled 
thrust chamber, the injector -inlet propellant enthalpy must include the heat 
transferred to the propellant from the chamber, and is not the propellant 
enthalpy at the engine inlet. 

It is desirable that the experimental value of a performance parameter 
be derived directly from the measured test data, without involving any 
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analytically based corrections whatsoever. Though this criterion is satisfied 
by the definition of the experimental vacuum specific impulse, it cannot be 
completely satisfied in the definition of experimental characteristic velocity 
because the basic definition of c* logically involves the stagnation pressure 
at the nozzle throat, where the mass flow is determined. Determination of 
the throat stagnation pressure to be used in the calculation of C' 1 * is compli- 
cated by the facts that (1) with multi -stream -tube flow, each stream tube has 
a different stagnation pressure, so that some sort of a mean value must be 
defined and used; and (2) it is not possible to measure directly the stagnation 
pressure of die flow in the throat of a rocket nozzle. 

The only pressure relevant to the stagnation pressure that can be mea- 
sured is the static pressure p s at some longitudinal station x on the wall 
of the thrust chamber. An analytically based correlation factor, (p^'Ps x )> 
then used to obtain a value for the mean effective stagnation pressure at the 
nozzle throat.^ 

For an ideal rocket engine with uniform perfect gas flow and a cylin- 
drical combustion chamber, the factor (pg/p Sx > is wel1 defined and P re " 
sented in the literature (see Fig. 8). This is the factor commonly used in 
computing characteristic velocity from the experimental data for all types of 
rocket motors. 

When the combustion chamber is not cylindrical and/or there is multi 
stream-tube flow, the ideal rocket engine factor does not apply. However, 
the DER computer program can be used for these cases to analyze the flow 
and to determine the static pressure-stagnation pressure relationships as 
well as the droplet evaporation and combustion efficiency. The results of 
a recent application of DER to determine the static pressure -stagnation 


^Static pressure measurements made in the injection region of the chamber 
may be in error by several percent due to local aspiration effects. A mid 
chamber or nozzle -entrance static pressure measurement is generally 
acceptable. 

14 The availability of the DER computer program opens up the possibility of 
defining and using a new mass -flow parameter based on a measured static 
pressure. The DER computer program would enable prediction of the 
static pressure at the location of the pressure measurement, and, thus, 
prediction of the analytical counterpart of this new parameter. However, 
this parameter would lack the generality of the present c*. 
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pressure relationship for multi -stream -tube flow with droplet-evaporation- 
limited combustion in a noncylindrical chamber are described in Ref. 4. 

An additional slight departure from the criterion of using only direct 
measurements in the calculation of experimental performance parameters 
occurs in the determination of the throat area to be used in calculating the 
characteristic velocity. Since it is not feasible to measure the throat area 
during motor operation, it is customary to measure this area under 'cold 
pre-test conditions and apply appropriate corrections for thermal expansion 
and stress, pressure, erosion, etc. These corrections are generally quite 

small. 

As a consequence of the considerations discussed above, it must be 

accepted that the characteristic velocity c* xp has an inherent uncertainty of 

2 to 3%. This arises largely from the uncertainty in inferring p^‘ from a 

measured static wall pressure p Sx> especially when there is multi-stream- 

tube flow through the nozzle, or when the combustion occurs in other than a 

cylindrical chamber. Despite this uncertainty in absolute level, c* xp can b< 

comoared with c* to give a useful first look at the magnitude of the 
pred 

losses involved in the thrust chamber combustion process. 

IX SIMPLIFIED FORM FOR ANALYTICAL PREDICTION 
OF CHARACTERISTIC VELOCITY 


A close analvtical approximation, following the JANNAF thrust 
chamber combustion and loss model, to the experimental characteristic 
velocity defined in the preceding section is 
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where c* , m* , and (O/F)? correspond to the local stream-tube 

i, vap i, vap L vap 

vaporized mass and mixture ratio at the throat. 
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Expression (26) is an adequate approximation to the overall mass -flow 
parameter as long as the liquid droplet fraction is a small part of the total 
mass flow. It is derived by summing the one -dimensional thermodynamic 
stream -tube throat areas and equating them to the total effective throat area 
derived by correcting the geometrical throat area for boundary -layer and 
two-dimensional flow effects. It is assumed that the individual stream-tube 
stagnation pressures can be replaced by a single mean stream -tube stagna- 
tion pressure, and that multi -stream -tube -flow sonic-point-displacement 
effects on the effective throat area are negligible. 

Expression (26) for predicted characteristic velocity can be recast as 
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pred 
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(27) 


The terms in the above expressions are defined and evaluated as follows 
where each factor r\ corresponds to a specific loss: 


(!) c ODE 

(0/F)avg 

is evaluated at the specified test conditions of chamber pressure, 
overall mixture ratio, propellant composition (including impu- 
rities), and inlet enthalpy. For a regeneratively cooled motor, 
the inlet enthalpy must include all of the heat transferred from 
the thrust chamber to the propellant. 

(2) The two-dimensional nature of the flow results in a curvature of 
the sonic surface and a net decrease in the effective one- 
dimensional flow area. This in turn results in a loss factor 
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Values of this factor are given as a function of the throat radius 
ratio in the summary chart prepared by Back and Cuffel (Figs. 9 
and 10 and Ref. 3). 

The kinetic loss arises from the fact that the combustion gas 
equilibrium cannot change fast enough to follow the changes in 
pressure and temperature as the gases expand through the nozzle, 
This effect is rarely important upstream of the nozzle throat, 
and thus has little effect on the characteristic velocity, the 
corresponding loss factor is 


E c f 


ODK 

< 0 /F >i*vat 


£ c * 


ODE 

(O/F)* ^ 

1 , vap 


-T) 

( — m? \ 

1. vap l 

r "t / 


c ODK 


(Q/F)avi 


c ODE 

( 0 /F)avgJ H 

inj 

size 


shape 


This factor must be calculated for the specific propellant, oper- 
ating conditions, and thrust chamber size and shape, as defined 
by Eq. (25), using the ODE and ODK computer programs. 

The vaporization loss factor accounts for the fact that the frac- 
tion of the propellant which has evaporated and reacted occupies 
most of the cross-sectional area of the nozzle throat. It also 
accounts for the fact that nonuniform droplet evaporation can 
cause vaporized mixture-ratio distributions (at the nozzle throat) 
which differ from the injected mixture-ratio distributions. This 
factor is defined as 
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The vaporized mass and mixture-ratio distributions needed to 

calculate n * are obtained from the DER computer program 
c vap 

calculation, which in turn requires as input a complete descrip- 
tion of the injected-mass, mixture -ratio, and droplet-size 
distributions. 

(5) The primary mixture-ratio loss arises from nonuniform injected 
mixture-ratio distribution. With the usual concave -downward 
variation of characteristic velocity with mixture ratio, the mass- 
averaged performance of a distributed mixture-ratio flow is less 
than the performance that would be achieved with the entire flow 
at its average mixture ratio. The factor accounting for this 
loss is 
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(6) The chamber injection-end heat-loss factor is defined as 
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This factor accounts for the convective heat lost to the system 
from the region (injector faceplate and a portion of the chamber 
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wall) upstream of the starting point for the boundary layer, and 
for radiation loss throughout the chamber; this is the heat lost 
in addition to that accounted for by the boundary-layer calcula- 
tions. This factor is less than unity only for the heat which is 
lost to the system or transferred to a regenerative coolant so 
that the measured injection enthalpy is increased; there is no 
net performance loss when the heat is picked up internally by the 
propellant (as in an injector faceplate and manifold) and returned 
to the combustion region (see Fig. 4). 

(7) The boundary -layer performance loss arises from friction and 
from cooling of the combustion products near the wall. The 
effect is usually to reduce the effective inviscid-flow area of the 
nozzle throat. The corresponding performance loss factor is 



The boundary -layer displacement thickness 6* is obtained 
from a boundary-layer program such as TBL, using TDK 
property output. 

The boundary-layer displacement thickness at the throat can become 
negative if the wall is highly cooled. Because of the accelerating flow through 
the throat, the boundary layer in this region is usually very thin, so that 
n , is usually very close to unity. 

c Sl 

Two of the above-defined losses, related to the uniformity of the mix- 
ing and the completeness of the evaporation, can be grouped together as an 
energy- release loss. Thus, 
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It is often convenient to consider n c g;R as a whole > sxnce lts individual 
components are related and interacting. 


The derivation and definition of the characteristic velocity efficiency 
factors are summarized in Table 4. 


X. CORRELATION OF EXPERIMENTALLY MEASURED AND 
PREDICTED CHARACTERISTIC VELOCITY AND 
COMPARISON WITH SPECIFIC IMPULSE 
PARAMETERS 


As is the case with specific impulse, the measured and predicted 
characteristic velocities can be compared to confirm or infer the magnitude 
of the individual losses. The losses can then be evaluated and a judgment 
made as to whether they are individually reasonable and acceptable, or 
whether thrust chamber or injector design changes should be made in an 
attempt to decrease the magnitude of those losses which are susceptible to 
control. The two-dimensional flow, kinetic, upstream chamber heat- 
transfer, and boundary-layer characteristic velocity losses are generally 
small. This leaves the energy-release losses (vaporization and mixture- 
ratio distribution) as the major factors in the overall characteristic velocity 
efficiency; these are also the losses most susceptible to individual control 
by changes in injector and chamber design. However, the control will often 
be at the expense of some other deliberately designed-for characteristic, 
such as low heat flux, chamber wall compatibility, low injector pressure 
drop, or small combustion-chamber size. 

A common application of characteristic velocity correlation is to 
obtain a first look at the energy-release efficiency, "backing it out of the 

experimental data and the more easily estimated losses, according to the 

15 

equation: 
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bracketed term in the denominator of Eq. (35) usually differs only 

slightly from unity, so that as a rough approximation = n c * ' 

ER TC 
exp 
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It must be remembered that vacuum specific impulse and the various 
specific impulse efficiencies are the ultimately important quantities in the 
design and operation of rocket motors. Characteristic velocity is a different 
performance parameter from specific impulse, and characteristic velocity 
efficiencies are not identical to the corresponding specific impulse effici- 
encies. However, the energy-release components of both the specific 
impulse efficiency and the characteristic velocity efficiency are related in 
that they are both ratios of mass-averaged values of their respective quanti- 
ties over the evaporated mass and mixture -ratio range compared to the 
value of the quantity at the average injected mixture ratio (see Eqs. 14 and 
30). It is clear that 


and 


sp 

ER 


'c* 

ER 


will both respond in the same manner to a deficiency in evaporated mass 
flow, and both will be influenced by the relative curvature of their respective 
values as functions of mixture ratio in the operating mixture -ratio region. 

In fact, if I S £ and c* were exactly similar functions of mixture ratio, 
the two energy-release efficiencies would be identical, but this condition is 
generally not exactly satisfied. 


Inspection of the relative curvature of the Isp vac anc * c * curves as 
functions of mixture ratio in the operating mixture -ratio region should give 
an indication of the relative response of the corresponding energy-release 
efficiencies to mixture -ratio distribution for a given propellant; it will 
usually be found that the specific impulse efficiency (at all nozzle area ratios) 
and the characteristic velocity efficiency are affected about equally by 


mixture -ratio distribution. Thus, 
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The difference between these quantities will usually be small. 
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The above relationships can be used in estimating 


n 


i 

sp 

ER 


from an experimentally determined value of q * • This is often a con- 


ER 


venient first step in a thrust chamber development or evaluation program, 


since the necessary testing can be conducted using a low -area -ratio thrust 
chamber in an ambient pressure environment and on a fixed test stand with- 
out thrust measurement. However, the slight extra effort of measuring 
thrust as well as propellant flow rate and chamber pressure, even in initial 
phases of testing and with low-exit-area-ratio thrust chambers, is generally 
worthwhile, as it enables determination of the ultimately needed specific 
impulse parameters. 


In conclusion, it should be stated that while characteristic velocity is 
a useful performance parameter in its own right (relating mass flow to 
throat area and chamber pressure), it is not a suostitute for specific 
impulse. Specific impulse data are an essential requirement for thrust 
chamber development and for performance correlation and prediction, and 
provision for obtaining the needed test data should be included in any such 
development program. 
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Table 1. Thrust chamber performance losses 


1. Chamber heat loss (upstream of boundary -layer start point) 
Z. Two-dimensional flow 

3. Kinetics: rate-limited equilibrium shift during expansion 

4. Vaporization: incomplete liquid droplet evaporation up to 

nozzle throat 

5. Mixture-ratio distribution at injection 

6. Boundary layer: friction and heat transfer at chamber and 
nozzle walls 


Table 2. Thrust chamber performance computer programs 


1 . 

ODE 

One "dimensional equilibrium combustion and 
nozzle flow 

2. 

DER 

Distributed energy . '-ease: mixture -ratio dis- 
tribution and dropler. evaporation up to nozzle 
throat 

3. 

ODK 

One -dimensional kinetic multi-zone nozzle flow 

4. 

TDK 

Two-dimensional kinetic multi-zone nozzle flow 

5. 

TBL 

Boundary layer, integral method 

6 . 

MABL 

Boundary layer, finite difference method, 
including mass addition 

BLIMP 
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Table 3. Development of approximate method for predicting vacuum 
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Fig 2. Internal processes in the real rocket thrust 
chamber 



Fig. 3. Experimental data needed for thrust chamber 
performance and evaluation 
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POINT WHERE H. . IS MEASURED) . 
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B. NO EFFECT. 


Effect on performance of heat losses from interior of 
thrust chamber 
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Fig. 7. Typical effect of throat wall radius of curva- 
ture on the two-dimensional specific impulse 
efficiency of conical nozzles 


JPL Technical Memorandum 33-548 



1.3 


A 

t 

A 

c 

p o 

p« 


: NOZZLE THROAT AREA 

CYLINDRICAL COMBUSTION CHAMBER AREA 

ISENTROPIC STAGNATION PRESSURE AT THROAT 

= STATIC PRESSURE AT INJECTION END OF CHAMBER 

= STATIC PRESSURE AT DOW NSTREAM END OF 
CYLINDRICAL COMBUSTION CHAMBER 
(COMBUSTION COMPLETE) 



Fig. 8. Effect of combustion chamber area ratio on 
chamber static pressure at injection end 
and downstream end of cylindrical combustion 
chamber 
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Fig. 9. Effect of nozzle throat radius of curvature on 
two-dimensional nozzle flow coefficient 
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APPENDIX A 

BOUNDARY-LAYER EFFECTS ON PERFORMANCE 

Boundary-layer effects in rocket nozzle flow arise from friction and 
from heat transfer at the nozzle surface. Consider a simple one- 
dimensional axisymmetric flow as shown in Fig. A-l. 

Far away from the wall, the flow is essentially inviscid and unaffected 
by the presence of the wall. In this region, it has "free-stream" values of 
pressure, temperature, velocity, and density. Near the wall the velocity 
and temperature profiles bend to match the values at the surface, as shown 
on Section AA and Section BB of Fig. A-l. The boundary-layer thickness 6 
is the distance from the wall at which there is no appreciable departure from 
free-stream conditions, either velocity or temperature; the static pressure 
is presumed to remain constant at the free-stream value throughout the 

boundary layer. 

The boundary-layer displacement thickness 6* is such that an inviscid 
free-stream flow extending from the axis out to the radial station (r - 6*) 
would have the same mass flow as the real flow, with boundary layer, extend 
ing from the axis to the wall. 

The boundary -layer momentum thickness 0 is a measure of the excess 
of momentum which the equivalent inviscid flow has over the real flow: an 
inviscid free-stream flow extending from the axis out to (r - 6* - 9) would 
have the same momentum flux as the real flow, with boundary layer, extend- 
ing from the axis to the wall. 

Then, the thrust of the real flow is given by 

(1) The thrust of the equivalent inviscid flow through a real nozzle 
of contour (r - 6*). 
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(2) Less the excess momentum flux of the equivalent inviscid flow. 

(3) Plus the force exerted by the exit pressure on the annulus of 
width 6* between the equivalent nozzle exit diameter and the 
real nozzle exit diameter. 

In algebraic form, taking account of the nozzle exit divergence angle, 
as shown on Fig. A*2, and neglecting second~order terms, the thrust of a 
real nozzle can be written: 

F* — p 

inviscid 
TDK 
€ ' 

where the subscript e' indicates that the calculations are made for inviscid 
flow through the equivalent nozzle contour. 

This can be written in the form below, where it is equivalent to Eq. (4) 
of the text: 
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Then, following Eqs. (5), (6), (7), and (8) of the text, the expression 

for A lap » to be used when the performance is referenced to an inviscid 
2 B L 

expansion through the geometrical nozzle contour, is developed. 

Thus, as in Eq. (7) of the text, 
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and from Eqs. (5) and (8) of the text 
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At this point, the boundary-layer specific in' pulse loss would be evalu- 
ated by the following steps: 

(1) Use the TDK computer program to determine the overall per- 
formance and the properties (p, u, p, T, c p , etc.) of the 
boundary flow stream tube. 

(2) Use the boundary stream-tube properties from TDK and a speci- 
fication of wall temperature distribution as input to a 
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boundary -layer computer program (TBL, MAbL, BLIMP) and^ 
determine the boundary-layer thicknesses (® e / r e ) and ( 5 */r g ). 

(3) Use the output of the previous two steps to determine the rates of 
change with nozzle area ratio of the quantities in the last two 
terms of Eq. (A-5); this can be done graphically. 

(4) Use the results of steps (1), (2), and (3), above, in Eq. (A-5) to 

calculate the value of A 2 I gp • 

BL 

The above calculation of the boundary -layer specific impulse correc- 
tion is essentially a "reference methodology" calculation, as it requires the 
use of both the TDK and one of the boundary -layer computer programs. It is 
the intent of this report to develop a "simplified" performance calculation 
procedure which avoids the use of long-running and therefore expensive -to- 
use computer programs like TDK and, if possible, the various boundary- 
layer programs. This objective has been achieved up to the point of 
computing the boundary-layer loss. (All of the other components and losses 
of the performance are calculated or approximated with the use of only the 
ODE, TDE, and ODK computer programs.) It is now necessary to attempt 
to find feasible and acceptable approximations and simplifications for the 
evaluation of the boundary -layer correction. 

As a first step in the simplification of Eq. (A-5), it is convenient to 
eliminate the mass flow rate in favor of the characteristic velocity. From 
Eqs. (26) and (27) of the text, the following equation can be written: 



rn 
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n c* r C &DE 

1Q ' (Q/F)avg 


(A-7) 


For further simplification, it can be assumed that the boundary-layer 
pressure, velocity, and density (and thus boundary-layer thicknesses) are 


16 


The combustion gas recovery temperature and heat-transfer coefficien , 
th wall temperature, and the heat flux into the wall must be an internally 
consistent set. It may take several iteration*.- with a heat-transfer pro- 
gram to achieve a heat balance and converge on the correct wall tempera- 
ture for a regeneratively cooled thrust chamber. 
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influenced much less by kinetic effects than by two-dimensional nozzle flow 
effects, so that TDE can be used instead of TDK to determine the fluid prop 
erties. 17 Additionally, the specific impulse in the last term of Eq. (A-5) 
can be expressed in terms of the one -dimensional equilibrium performance 
and the applicable efficiency factors, and the efficiency factors which are 
insensitive to nozzle area-ratio variations can be removed from the partial 
differentiation operation. With these substitutions and approximations, 

Eq. (A-5) for the boundary -laye r specific impulse correction becomes 
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As expressed by Eq. (A-8), the boundary-layer specific impulse correction 
would be computed by the following procedure: 


(1) Use ODE to computer I and c*. 

sp 

(2) Use ODK to determine ni gp . 

kin 


These operations would 
already have been 
performed. 


^ 7 A two-dimensional solution is required because the nozzle wall curvature 
can significantly affect the pressure distribution along the boundary layer, 
and this enters into the evaluation of the boundary-layer displacement and 
momentum thicknesses. A TDE option is being prepared as part of a 
revision of the TDK computer program. 
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( 3 ) 


Use TDE to compute boundary-layer stream properties (p e , P e> 
V T e’ C p’ CtC * ** 

(4) Use a boundary-layer computer program to compute 0 g and 6*. 

(5) Determine ’ll from available design charts and calculate or 

2D 

estimate 


. n T . and r\ „„ 

I 9 I c nrr 

sp sp 


HL 


ER 


(6) Determine graphically the slopes occurring in the second and 
third terms of Eq. (A-8). 

Examination of the relative magnitude ot the three terms of Eq. (A 8) 
shows that the first term is predominant. The second and third terms are of 
opposite sign and become relatively negligible at high nozzle area ratio. 

Also, at high nozzle area ratio. the relative effect of the (b*/Q Q ) factor in 
the first term decreases. At low area ratio, corresponding to "sea level" 
testing, all terms in Eq. (A-8) will have to be included in the calculations. 

In the numerical evaluation of the first and second terms of Eq. (A 8,, 
it is helpful to recognize that 

(1) The term U (p /p*) (p^ u^/p )] is relatively invariant, having a 
value of about 1. 25 at e = 2 and increasing to about 1 . 8 at € - 100 
(see Fig. A-3). 

(2) The term (p u^/p g ) varies from y at the nozzle throat to from 
20 to 50 at « = 100, depending on the value of y (see Fig. A-4). 

The evaluation cf the boundary-layer specific impulse loss from 
Eq. (A-8), as described above, is a simplified direct calculation procedure 
which uses the most economical of the computer programs giving the needed 
output information, and is still based on the actual combustion product prop- 
erties and the actual thrust chamber and nozzle size and shape. 
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An approximate boundary-layer specific impulse loss evaluation 
method is described in Ref. 1, Appendix B. This method uses charts based 
on calculations using the TBL boundary-layer program for typical thrust 
chamber nozzle shapes and over a range of wall temperatures and combus- 
tion gas properties. This method can be used to obtain rough estimates of 
the boundary-layer loss. 
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Fig. A- 1 . Boundary-layer parameter definitions 



Fig. A- 2. Effect of boundary layer on nozzle geometry 
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Fie A- 3 Perfect gas values of the thermodynamic function U (P e /Po )( p e U e^ P e ^ used 
' latine performance loss due to boundary-layer effects 






APPENDIX B 


FLOX-MMH PROPELLANT DATA 

This appendix is presented as an example of the precalculations that 
can be performed for a given propellant, using the ODE and ODK computer 
programs. The propellant chosen as an example is the Flox (88% F.,-12% O-,) 
oxidizer with MMH fuel. The enthalpy of the propellant components (injection 

1 8 

enthalpy) is as follows: 


Propellant 

component 

Molecular 

weight, 

g/g-mol 

Enthalpy, 

cal/g-mol 

Temperature, 

K 

Specific 
heat, 
cal ; g-K 

MMH (CH 6 N 2 ) 

46. 08 

12700 

298. 15 

0. 6928 

f 2 

38 

-3100 

90. 20 

0. 363 

°2 

32 

-3080 

85. 20 

0. 450 


The performance calculations were made for a chamber pressure of 
100 psia. Kinetic effects were determined at the scale of a 600-lbf thrust 
chamber (nozzle throat diameter = 2. 00 in. ). 

The ODE computer program is used to determine enthalpy of the mixed 
propellant, the combustion temperature, the characteristic velocity, and the 
vacuum specific impulse as a function of propellant mixture ratio. The 
results of these calculations are shown in Figs. B-l and B-2. 


18 A set of recommended values of enthalpy, specific heat and density for 
various propellants is given in Appendix 2 of the Minutes of the Sixth Meet 
ing of the JANNAF Performance Standardization Working Group, 

October 26-27, 1972, and is reproduced here as Table B-l. 
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Equivalent calculations are made with the ODK computer program, and 
the ratio of ODK to ODE values for characteristic velocity and specific 
impulse is plotted as kinetic efficiency in Fig. B-3. 

The ODE calculations were repeated with higher and lower values of 
the propellant enthalpy in order to determine the change in characteristic 
velocity and specific impulse due to cnanges in injection enthalpy. The 
influence coefficients obtained from these calculations are given in Fig. B 4. 
These influence coefficients can be used, as described in the text, in account - 
ing for differences in the total injection enthalpy from the values shown in the 
above tabulation, and in accounting for heat loss to the injector and chamber 
upstream of the boundary-layer attachment point. 

In the absence of a good set of boundary -layer calculations for this 
propellant in the thrust chamber configuration shown in Appendix C, the 
boundary-layer specific impulse efficiency decrement was estimated as 
shown in Fig. B-5 for use in the sample performance data analysis. 
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REPRODUCIBILITY OF THE ORIGINAL PAGE IS POOR. 


Table B-l. Recommended heats of formation for propellants 

______ _____ Assign* Enthalpy Temperature Density 


Empirical Formula 


Compound 

Chlorine Pfhtafluorlde (L) 

l inni _ 0*\ 


Cp(cal/gm-°K) • 0. 1 941 ♦ 0.4824 * 10‘ 3 {T,°K) 
over ra^ge 223-323*K (Ref. 2) 

Chlorine Trifluoride (l) 

Cp(BTU/lb* t> K) ■ 0.30939 - 1.672 x 10" 4 (T,°R) ♦ 3.0870 x 10 *(T, 9 R)* 
Over range 355-500°R (Ref. 5) 

Dlborane (L) 

Fluorine (1) 

Cp(c#l /mol -*K) * 13.8 • 82°K (Ref. 2' 

Hydrazine (t) 

Cp(cal -gm^K) - 0.7356 9 298°* , _ 6 2 

Cp(cal/gm-°K) • 0.88415 - 1.3949 x 10 3 (T.*K) ♦ 3.0074 * 10 (T, R) 

(Ref. 3 anhydrous hydrazine) 

2 Hydrogen (1) 

0^ Hydrogen Peroxide (L) 

Cp(cal/gm-*K) • 0.628 average between 273 and 300"K 

0. .... 901 Hydrogen Peroxide (l) 

2 1 . 8266 — - 


101 H 2 0 (1) 


Cp(cal/^n*°K) • 0.660 average between 273 and 300°K. empirical 
formula assumes one mole solution 

;h 4 Methane (L) 

>« w Monomethyl hydrailne (L) 

o 2 

Ref. 2 recommends aH° 29g • 12.7 keel /mol 
Ref. 3 recommends aH* 29 ^ • 13.106 keel /mol 

Cp(cal/**-*K) ■ 0.6528 - 1.7284 x 'l0' 5 (T.*K) ♦ 3.9142 x 10* (T.°R) 
over range 226-360*R (Ref. 3) 

Nitric Acid (l) 

- * 090 . .o; s (T.;^ 

no temperature range given for heat capacity equation (Ref. 2) 

V»W N 0.9365°2.*» ™ * C,d " FHA ( ° 

|4t HN0j/14l NOj/21 MjO. empirical formula asstmes one mole solution 
n 2 0 4 Nitrogen Tetroxide (L) 

Cp(ITU/lb- # R) • 0.3698 • 536. 7*R 4 

Cp(»TU/1b-**) • 0.2*668 ♦ 2.2880 x 10 *(T. B R) Rtf. 5 
from 477 to 524*R 

0 2 Oxygen ( L ) 

Cp(ce1/gm-*K) - 0.4050 • 90.18^ n.b.p. (Ref. 2) 

0 F ? Oxygen 01 fluoride (l) 

C ? h # n 2 unsym Olmtthylhydrellne (L) UCAM 

Ref. 3 recommends aHf^ • 12.339 kctl/mol 
Cp(ca1/*n-*R) • 0.653 9 298*il (Ref. 2) 

Cp(ce1/9*-*K) • 0.4071 ♦ 8.838 x 10’ 4 (T,*K) 

Over range 216-335 R by Re*. 3 but values considered provisional 

WYjniYo Mi 501 ’ J ™* 1 V« a) 

Enthalpy obtained by molar addition of N^i (0.6522 mole) end 
UOAM (0.3478 mole) plus heat of mixing 
Cp(cal/^*R) • 0.732 t m-R 4 


Temperature Density References 

“ic g/cc Enthalpy Density 


-60.5*6.0 

-46.4 

298.15 

1.779 

-44.45*0.8 

-48.08 

298.15 

1.807 

♦4.97*4.0 

♦ 17.94 

180.59 

0.4371 

-3.098*0.02 

-8.153 

85.02 

1 .505 

♦12.05 

♦37.60 

298.15 

1.0037 

-2.154*0.02 

-106.8 

20.27 

0.0709 

-44 . 88*0.02 

-131.9 

298.15 

1.44 

-45.01 

-144.07 

298.15 

1.4136 

-21.39*0.10 

-133.3 

111.66 

0.4239 

*12.9 

*28.0 

298.15 

0.8702 

-41.46*0 10 

-65.79 

298.15 

1.5027 

-36 48 

-63.81 

298.15 

1.55 

-4.68*0.4 

-5.08 

298.15 

1.431 


-3.102*0.02 

-9.694 

90.18 

1.149 

•8.38*0.90 

-15.58 

127.8 

1.521 

♦11.9 

♦ 19.8 

298.15 

0.7861 

♦12.31 

♦29.45 

298 15 

0.8987 


CP cel/p*R) • 0.5124 * 7.3824 x 10‘ 4 <T.*I) ever range 27»-123*l 
(Ref. 3), empirical fer»**1e as**** ene mole solution 
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MIXTURE RATIO (O Ft 


Fig. B 


1. ODE reference performance for Flox (88-12)-MMH (enthalpy, 
chamber temperature, characteristic velocity) 
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MIXTURE RATIO (O/F) 


Fig. B-Z. ODE reference performance for Flox (88-12)-MMH (vacuum 
specific impulse vs nozzle area ratio) 
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AREA RATIO t 


Fig. B 


5. Boundary-layer loss estimate for 600-lbf thrust 
chamber 
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APPENDIX C 


EXAMPLE CORRELATION AND EXTRAPOLATION 
OF FLOX-MMH TEST DATA 

The methods of thrust chamber performance correlation and 
prediction developed and discussed in the preceding portions of this report 
are applied here, as an example, to experimental data obtained at the Jet 
Propulsion Laboratory. 

The experimental data was obtained during tests leading to the develop- 
ment of a 600 -lbf vacuum thrust rocket motor utilizing Flox (88-12)-MMH as 
propellant. The test motor had a copper heat-sink thrust chamber with a 
nozzle exit-area ratio of « = 2. 50, and was operated at a chamber pressure 
(nozzle throat stagnation pressure) of approximately 100 psia. 

Data from a number of tests made over a range of mixture ratio were 
plotted versus mixture ratio, and the performance values used here for the 
example correlation were read at a mixture ratio of (O/F) - 2. 20 from 
curves faired through the plotted data points; these performance values are 
given in the labeled box in Table C-l. 

Table C-l shows performance data for two different thrust chamber 
lengths, corresponding to L* values of 18 and 40 in. Comparison of these 
two sets of data and the corresponding performance efficiency factors 
clearly shows that the combustion and energy release is more complete in 
the longer thrust chamber. The right-hand column of Table C-l shows the 
extrapolation of the efficiency factors obtained by correlation of the € - 2. 5 
test data to a thrust chamber with a nozzle exit-area ratio of c - 60, and the 
predicted specific impulse at this area ratio. 
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The correlation of the experimental data and the prediction of the 
performance with the * = 60 nozzle are obtained by the following steps: 

(1) Correlation of € = 2. 5 test data 

(a) Enter in the table the "reference" equilibrium performance 
of the propellant at the over?\l average mixture ratio, as 
given in Appendix B, Figs. B-l and B-2. 

(b) Correct this "reference" performance for the difference of 
the actual propellant component injection temperatures and 
enthalpies, compared to the "reference" values, as listed 
in Appendix B. This is done using the injection enthalpy 
influence coefficient curve (Fig. B-4). The sample calcu- 
lations are given as Appendix D. 

(c) Compare the measured delivered performance with the ODE 
equilibrium performance for the propellant at the actual 
injection enthalpy to obtain the thrust chamber performance 
efficiencies (Eqs. 2 and 25 of the text). 

(d) Compute the performance efficiency due to the heat loss 
from the injector face and the portion of the chamber 
upstream of the boundary -layer attachment point, using 
Eqs. (10) and (32) of the text and the influence coefficients 
from Fig. B— 4. In this case it is assumed that 3 in. of 
chamber length is involved, and that the heat flux in this 
region is q = 1.0 Btu/in?-sec. The details of this calcu- 
lation are given in Appendix E. 

(e) The efficiency factor for two-dimensional flow through the 
nozzle is read directly from Fig. 5 or 6 for specific 
impulse, and obtained from Fig. 9 and Eq. (28) of the text 
for characteristic velocity. 

(f ) The efficiency factors resulting from kinetic loss at the 
overall average mixture ratio are read directly from 
Fig. B-3. 

(g) The effect of friction and heat loss from the attached 
boundary-layer region on specific impulse is taken from 
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4 * 


the estimated curves of Fig. B-5. The effect of the 

boundary -layer loss on characteristic velocity is assumed 

to be = 1. 002, corresponding to the assumption that 

(6* /r ) = 0. 001 (cf Eq. 33 of the text), 
t t 

(h) The energy-release efficiency is then "backe^ out" from 
the measured overall thrust chamber efficiencies and the 
other component efficiencies evaluated above, using 
Eqs. (23) and (35) of the text. 

(2) Extrapolation to performance at e = 60 

(a) The predicted specific impulse at t = 60 for the propellant 
at its actual injected enthalpy is found as before, using 
Figs. B-l, B-2, and B-4. The associated calculations are 
given in Appendix D. 

(b) The injector and chamber upstream-end heat loss specific 
impulse performance efficiency is obtained as before from 
Eq. (10) and the influence coefficient given in Fig. B-4. 
The calculational details are given in Appendix E. 

(c) The two-dimensional flow effect on specific impulse at 

€ = 60 is obtained directly from Figs. 5 and 6 of the text. 

(d) The kinetic efficiency factor for specific impulse at « = 60 
is obtained directly from Fig. B-3. 

(e) The boundary-layer effect on specific impulse at e = 60 is 
obtained from the curve in Fig. B-5. 

(f ) The energy release efficiency at £ = 60 is conservatively 
assumed to be the same as at t - 2. 5. 

(g) The overall thrust chamber specific impulse efficiency at 
t = 60 is then calculated from the component efficiencies, 
using Eq. ( 1 6) of the text. 

(h) The predicted, or extrapolated, vacuum specific impulse 
at £ = 60 is obtained by multiplying the ODE specific 
impulse by the overall thrust chamber specific impulse 
efficiency factor calculated above, as per Eq. (15) of the 
text. 
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The results of the data correlation are summarized in the first two 
columns of Table C-l; it is found here that the characteristic velocity energy- 
release efficiency is very nearly equal to the overall thrust chamber charac- 
teristic velocity efficiency. For this particular example, too, the specific 
impulse energy release efficiency at (O/F) = Z. 2 and at e - 2, 5 is very close 
to the characteristic velocity energy-release efficiency. 


74 


JPL Technical Memorandum 33-548 


Table C-l. Test data and efficiency factor correlation for tests of 600-lbf 
Flox (88-1 Z)-MMH thrust chamber with e = 2. 5 nozzle area ratio 
and extrapolation to performance at £ - 60 nozzle area ratio 



Derived : backed out 




APPENDIX D 


CORRECTION OF ODE PERFORMANCE FOR CHANGE 
IN PROPELLANT INJECTION TEMPERATURE 

The change in performance due to changes in propellant injection 
temperature is obtained by determining the equivalent change in propellant 
injection enthalpy, and then using the precalculated enthalpy change influence 
coefficients given in Fig. B-4. 

The "reference” ODE performance was computed for the propellant 
component temperatures given in Appendix B. These temperatures and the 
average temperatures of the propellant during the test program are given 
below: 


Propellant 

component 

Standard 

temperature, 

K 

Test 

temperature, 

K 

AT, K 

MMH (CH 6 N* 2 ) 

298. 15 

308. 0 

+10.0 

F 2 

90. 20 

82. 0 

1 

00 

°2 

85. 20 

82. 0 

+3. 2 


The total change in enthalpy, per unit mass of total propellant, is 
obtained by summing mass fractions of the propellant components times their 
specific heat and temperature change. Thus, 



for each propellant component. 
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The mass fractions of the propellant components at a mixture ratio of 
(O/F) = 2. 20 are obtained as follows: 

(O/F) = r 


oxidizer 


— - — = 4^4 = 0 . 688 

(r + 1) 3.2 


Since the oxidizer is 88% F^ and 12% 


= 0 . 88 x 0. 688 = 0. 605 


(a 


= 0. 12 x 0. 688 = 0. 0825 


\ rot 


(r + 1) 


= 0.312 


■ft) 


MMH 


Then, using the specific heats given in Appendix B, 


AH = AH p + AH Q ^ + AH mmh 


= [0. 605 x 0. 363 x (-7. 8) + 0. 0825 x 0. 450 x 3. 2 
+ 0. 313 x 0. 6928 x 10. 0] cal/g 

= + 0. 572 cal/g of total propellant 
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Then, from Fig. B -4 for (O/ F) - 2. 2, 



1. 25 x 10 
1. 33 x 10 


1 

cal/ g 

1 

cal/g 


1. 765 x 10 


1 

cal/g 


and the corresponding values of performance change are 

= 1.25 x I0* 4 x 0.572 = 0.0000715 

c* 

/ AXg P ) = 1. 33 x 10" 4 x 0. 572 = 0.0000760 

' Is P 4 = 2.5 


/ AIs p \ = 1. 765 x 10’ 4 X 0. 572 = 0.000101 

' Is P/ € = 60 

These values are completely negligible, so Table C-l shows the ODE 
performance corresponding to the injection enthalpy as being identical to the 
performance at the reference enthalpy. 
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APPENDIX E 


CALCULATION OF INJECTOR REGION HEAT LOSS EFFECTS 

Heat loss to the downstream portion of the combustion chamber and in 
the contraction and expansion portions of the nozzle is accounted for m the 
boundary-layer performance efficiency loss factor. There is an effective 
starting point, or boundary -layer attachment point, beyond which the 
boundary -layer processes will adequately predict the measured heat loss. 
Upstream of this effective starting point, in the injection region, there is 
considerable large-scale turbulence, and the heat transfer is not described 
by the normal boundary -layer relationships. The effect on thrust chamber 
performance of the heat loss from the region upstream of the effective 
attachment point of the boundary layer must be accounted for separately 
from the boundary-layer loss accounting. 

This appendix illustrates the method of calculating the effect of 
injector-region heat loss on thrust chamber performance. The configuration 
of the test chamber is shown in Fig. E-l. It is assumed that the effective 
attachment point of the boundary layer is 3 in. downstream of the injector 
face, and that the injector face and the chamber wall in the region upstream 
of the attachment point have a heat flux of 1. 0 Btu/in. -sec. 

The total surface area included for injection-region heat loss is 

A = — X (3) 2 + tt x 3 - 3 = 7. 075 + 28. 3 = 35 in? 

4 

The heat flux in this region is 

Q = 35 in? • 1 . 0 Btu/in?-sec = 35 Btu/sec 
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The total propellant mass flow into the chamber is given by 


m. - 


p Q A 100 (lbf/in?) X j (2) 2 (in?) X 32. 174 (lbm-ft/lbf-sec 2 ) 


t n„* c* 


c| c ODE 


0. 95 X 6875 (ft/sec) 


= 1.55 lbm/sec 


where the value of r) ^ = 0. 95 is assumed. 

C TC 

It is assumed that, because of the turbulence in the injector region, the 
heat is lost uniformly from the total propellant. Then the equivalent enthalpy 
change of the propellant is 


AH 


HL 


= A = 35 Btu/.ec. = 22 . 6 Btu/lbm 

1.55 lbm/sec 


m* 


= 12. 55 cal/ g 


From Fig. B-4 the injection enthalpy change performance influence 
coefficients at (O/F) = 2. 2 are 


1 Ac* \ _ , 

c* AH j ~ 


25 X 10 


-4 1 


:al/g 


i AI 

_i S£ 

I AH 
sp ! € =2. 5 


>- 33 * 10 ' 4 ^u 


g 


i AI 

1 sp 
I "AH 
sp / e =60 


>- 765 x 10 " 4 ^T7i 


The corresponding performance changes are 


= 1.25 x 10 -4 x 12.55 = 0.00157 


r AI 


sp 


= 1. 33 x 10 -4 x 12. 55 = 0.00167 


t =2. 5 
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AI \ .4 

®£j = 1.765 X 10 4 X 12. 55 = 0.00221 

W««60 


and the corresponding performance efficiency factors are 


V* 


HL 


sp 

HL 


(1 - = 0.9984 


/ AI \ 

= ( 1 - - S P ) = 0. 9983 at € = 2.5 
\ sp ' 


sp 

HL 


■( 


AI v 

1 ££ \ = 0.9978 at = 60 


I 

sp / 
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NOMENCLATURE 


AJ« 

l 


HP 

A t 

C 

C D 

2D 

* amb 

F 

vac 
AF 


BL 


AF 


drops 


H. . 
inj 


sp 

vac 


m?' 

1 


m:> 


i, drop 


m* 


l, vap 

(0/F)avg 

(O/F). . . 

i. inj 

(O/F).* 

i, vap 


P(S. 


geometric nozzle -exit area under firing conditions 
throat area of individual stream tube 
nozzle -exit lip area 

geometric nozzle throat area under firing conditions 
characteristic velocity — mass-flow parameter 
discharge coefficient for throat curvature effects 

thrust measured with external pressure p a 
thrust in vacuum environment 

boundary-layer correction to calculated thrust 
droplet contribution to total thrust 
injection enthalpy of propellant 
vacuum specific impulse 

mass-flow rate in one stream tube 

mass -flow rate of unevaporated droplets remaining at 
stream-tube throat 

evaporated mass -flow rate at stream -tube throat 
total mass -flow rate 

overall average injected mixture ratio of propellant 
injected mixture ratio in one stream tube 

mixture ratio of evaporated propellant at stream -tube throat 
ambient pressure 

pressure in boundary layer at nozzle exit 

static pressure at wall of combustion chamber at axial 
location x 

stagnation pressure at throat of stream tube (isentropic) 
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NOMENCLATURE (contd) 


Q 




n c * 

kin 


n r * 


(0/F)dist 


c ^ 
C TC 

exp 


n r * 


TC 

pred 



average stagnation pressure at throats of stream tubes 
(defined) 

rate of heat transferred to chamber upstream of boundary- 
layer attachment point (see subscript HL) 

radius of curvature of nozzle throat 

geometric radius of nozzle exit 

geometric radius of nozzle throat 

velocity, stream property at nozzle exit near wall 
divergence angle of nozzle wall at exit 
ratio of specific heats of gas 

displacement thickness of boundary layer at nozzle exit 
displacement thickness of boundary layer at throat 

correlation coefficient for experimental vs. predicted 

performance 

nozzle exit-area ratio 

characteristic velocity efficiency accounting for boundary- 
layer friction and heat-transfer effects 

kinetic characteristic velocity efficiency 

mixture-ratio distribution characteristic velocity efficiency 


C ODE 

(0/F)avg 


c* 

pred 

<0/F)avg 


characteristic velocity efficiency accounting for heat lost to 
the system upstream from the boundary “layer attachment 
point 


84 


JPL Technical Memorandum 33-548 


NOMENCLATURE (contd) 


n 

n 


C * 

vap 


sp 

kin 

n I 

sp 

(0/F)dist 



sp 

TC 


exp 


sp 

TC 

pred 


sp 

HL 


sp 

vap 


^1 


sp 

2D 


An, 


sp 

BL 


droplet vaporization characteristic velocity efficiency 
two-dimensional characteristic velocity efficiency 
kinetic specific impulse efficiency 


mixture —ratio distribution specific impulse efficiency 


I 

sp 

vac 

exp 

I 

sp 

ODE 

(0/F)avg 


I 

sp 

vac 

pred 

_ 

sp 

ODE 

(0/F)avg 

specific impulse efficiency accounting for heat lost to the 
system upstream from the boundary -layer attachment point 

droplet vaporization specific impulse efficiency 


two -dimensional -flow specific impulse efficiency 


specific impulse efficiency increment due to boundary-layer 
friction and heat transfer 

momentum deficiency thickness of boundary layer at nozzle 
exit 

density, stream property at nozzle exit near wall 
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NOMENCLATURE (contd) 


Subsc ripts 
BL 
drops 
e 

exp 

ER 

i 

inj 

kin 

ODE 

ODK 

(O/ F )avg 
(0/F)dist 
pred 
HL 

rad 

TDK 

vap 

x 

€ 

ZD 


boundary -layer loss effects 
unevaporated liquid droplets 
nozzle exit 

experimental value, based on measured data 
effects due to incomplete energy release 
individual stream-tube values 

corresponding to injected mass or mixture ratio 
effects due to finite reaction rates 

calculated v^ing one “dimensional equilibrium computer 
program 

calculated using one -dimensional kinetic computer program 
corresponding to overall average mixture ratio 
effects due to mixture-ratio distribution 
analytical prediction of experimental value 

effects due to heat lost to injector face and to chamber wall 
upstream of boundary-layer attachment point 

radiation effect 

calculated using two-dimensional kinetic computer program 
corresponding to local evaporated mass or mixture ratio 
axial position in thrust chamber 
corresponding to the nozzle exit-area ratio 

corresponding to the equivalent inviscid-flow nozzle area 
ratio 

two-dimensional flow effects 
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